Error Threshold for Color Codes and Random 3-Body Ising Models 
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We study the error threshold of color codes, a class of topological quantum codes that allow a direct im- 
plementation of quantum Clifford gates suitable for entanglement distillation, teleportation and fault-tolerant 
quantum computation. We map the error-correction process onto a statistical mechanical random 3-body Ising 
model and study its phase diagram via Monte Carlo simulations. The obtained error threshold of pc — 0.109(2) 
is very close to that of Kitaev's toric code, showing that enhanced computational capabilities do not necessarily 
imply lower resistance to noise. 

PACS numbers: 03.67 .Lx, 75.40.Mg, 03.67.Pp, 75.50.Lk 



Protecting quantum states from external noise and errors is 
central for the future of quantum information technology. Be- 
cause interaction with the environment is unavoidable, active 
quantum error-correction techniques based on quantum codes 
have been devised to restore the damaged quantum states from 
errors caused by decoherence \ lkQi\. These approaches are, in 
general, cumbersome and require many additional quantum 
bits, thus making the system more error prone. An imagina- 
tive and fruitful approach to quantum protection is to exploit 
topological properties of a system, e.g., by using the nontriv- 
ial topology of a surface to encode quantum states at the logi- 
cal level yl]. Topology is thus considered as a resource, much 
like entanglement is a resource for quantum information tasks. 
Topological quantum computation is the combination of these 
two resources with the aim of winning the battle against de- 
coherence. These topological quantum error-correcting codes 
are instances of stabilizer quantum codes jill, in which errors 
are diagnosed by measuring certain check operators or stabi- 
lizers. In topological codes these check operators are local, 
which, in practice, is an important advantage. Moreover, er- 
ror correction has a deep connection to random spin models 
in statistical mechanics and lattice gauge theories l^. 

One of the original motivations for introducing surface 
codes was to achieve error protection at the physical level 
through energy barriers that would remove the need for exter- 
nal recovery actions. Only the application of strong magnetic 
fields (compared to the topological coupling) destabilizes the 
topological phase However, several studies fHSSSESl 
and a rigorous proof jH] have shown that the toric code (TC) 
is not stable against thermal excitations, except in four dimen- 
sions Si]. 

Therefore, the study of active error correction in topolog- 
ical codes [511 is fully justified. Ultimately, the goal is not 
only to achieve good quantum memories, but also to per- 
form quantum computations with them. In this regard, the 
TC fstl is somehow limited since it only allows for a conve- 
nient (transversal) implementation of a limited set of quantum 
gates: Pauli gates of X and Z type, and the CNOT gate. To 
overcome this limitation, topological color codes (TCC) have 
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FIG. I: (Color online) p-Tc phase diagram for the random 3-body 
Ising model. For p > Pc ~ 0.109 the ferromagnetic order is lost. 
The dotted line is a guide to the eye, the black circle represents the 
analytically known transition temperature of the 2D Ising model. The 
blue (solid) line represents the N-line. In the regime marked by a 
dashed line the exact determination of Tc(p) is difficult. 



been introduced 111 IL 11211 . Using TCC, it is possible to im- 



plement the whole Clifford group of quantum gates and thus 
realize quantum distillation, teleportation, etc. Notice that, al- 
though we use the mapping of Ref. |5], there is a difference 
regarding the issue of types of homology involved: our model 
has a colored homology, while the Kitaev model has a sim- 
ple homology. As a result, unlike the standard Ising model, 
the resulting statistical mechanical model has 3-body interac- 
tions with a value of pc a priori unknown, thus motivating the 
present study. 

The question arises as to whether the wider computational 
capabilities of TCCs imply a lower resistance to noise. We 
address this problem and show that the (error) threshold value 
is P c = 0.109(2), which is comparable to values for the TC 



|13|,[14|,[15ll. To compute we derive a statistical mechanical 



model describing the error-correction process; a random 3- 
body Ising model, with (classical) spins located at the vertices 
of a triangular lattice. In addition to thermal fluctuations, the 
mapping requires the introduction of quenched randomness 
to the sign of the interactions that correspond to faulty bits. 
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FIG. 2: (Color online) Lattice for the TCCs with 3-colored vertices. 
Physical qubits of the error-correcting code correspond to triangles 
(stars mark the 'boundaries' of the sets of triangles displayed), (a) 
Boundary of a vertex v. The stabilizer operators , have sup- 
port on the corresponding qubits. (b) Error pattern in the form of a 
string net. The three vertices that form its boundary are all the in- 
formation we have to correct the error, (c) Two error patterns with 
the same boundary. Because together they form the boundary of the 
three vertices marked with a circle, they are equivalent. 



One can then study the p-T^ phase diagram of the model, see 
Fig. [H where p is the probabihty for wrong-sign coupUngs to 
appear. For low T, p the model orders, which corresponds to 
feasible error correction. The critical for error correction is 
recovered from the critical p along the Nishimori (N) line III6II 
in the p-T plane. 

The disordered 3-body Ising model on a triangular lattice 
has not been studied before. However, in the absence of ran- 
domness it is known to have a different universality class than 
the standard Ising model, but with the same critical temper- 
ature lUrtl. Furthermore, the critical exponents can be com- 
puted exactly {v = a = 2/3), which allows us to test the 
numerical results in the p = limit. 

Topological color codes. — To construct a TCC C we start 
from any two-dimensional (2D) lattice in which all plaquettes 
are triangles and vertices are 3-colorable, such that no link 
connects vertices of the same color. The lattice is embedded 
in a compact surface of arbitrary topology. Since informa- 
tion is encoded in topological degrees of freedom, the code is 
nontrivial only when the topology of the surface is nontrivial, 
e.g., a torus of genus g > 1. So far, color codes have been 



introduced in the dual lattice (2-colex 111 111 ). Here we prefer to 
work in the triangular lattice to have a more direct mapping, 
see Fig. [2] 

We consider a physical system with a qubit at each lat- 
tice triangle, and introduce the following vertex operators 
that generate the stabilizer group of C. For each vertex v 
we have two types of operators which correspond to Pauli 
operators of X or Z type, i.e., Xy := 



A:veA 



Xa and 



Zv '■= ^augA^a- Thus, a vertex operator acts on all 
nearby triangles, see Fig. [2] Vertex operators pairwise com- 
mute and square to identity. The code C is defined as the sub- 
space with Xy = Zy = 1 \/v. perform error correction 
one measures vertex operators. The resulting collection of ±1 



eigenvalues is the error syndrome. 

Error correction. — Color codes have a structure with sta- 
bilizer generators which are either products of X or Z Pauli 
operators, but not both. This allows us to treat bit-flip and 
phase errors separately, making the procedure classical: X- 
type (Z-type) errors produce violations of Z-type (X-type) 
vertex operators. Without loss of generality, let us consider 
the bit-flip case, that is, errors of the form X^; (S* asb -^^^ 
where E is the subset of triangles that suffered a bit-flip. Let 
dE be the collection of vertices that are part of an odd num- 
ber of triangles in E, i.e., the boundary of a set of triangles E 
is chosen so that the error Xe gives rise to a syndrome with 
Zy = — 1 at those vertices v S dE, see Fig.|2] In trying to 
correct the error, we apply to the system bit-flips Xe' with the 
same boundary, dE' = dE. This is only successful as long as 
Xe'Xe =■ Xe+e' is an element of the stabilizer group. Ge- 
ometrically, D = E+E' is a cycle: its boundary dD is empty. 
Given a vertex v, let dv be the subset of triangles meeting at v. 
We say that _D is a boundary if D = J2v ^'^^ some subset 
of triangles V. In that case, Xd is an element of the stabilizer 
group. Thus, error correction is successful whenever D is a 
boundary, i.e., if D has trivial homology. In that case the real 
error E and the guessed error E' belong to the same homology 
class. 

Mapping to a random 3-body Ising model. — We consider 
a standard error model based on stochastic errors in which 
phase errors Z and qubit bit-flip errors X are uncorrected 
and occur with probability p at each qubit. We focus on the 
correction of bit-flip errors. 

Let P{E) be the probability for a given set of bit-flip er- 
rors E. Up to a p-dependent factor, P{E) oc [p/ (1 — p)]''^'. 
We may also consider the total probability for the correspond- 
ing homology class E of errors, P{E) P{E + D), 
where D runs over all boundaries. If we measure a syndrome 
dE, then the probability that it was caused by an error in the 
homology class E is 



P[E\dE) 



P{E) 



(1) 



where the Di are representatives of the homology classes of 
cycles jslj. Then, error correction is achievable if in the limit 
of infinite system size we have J2e ^i^) P{E\dE) 1. 
That is, p < pc if for those syndromes which have a nonnegli- 
gible probability to appear the error can be guessed with total 
confidence. 

Following Ref. (E], we set cxp(— 2i^) := p/ (1 — p) (with 
K = J/T, T the temperature) for the N-line so that P{E) oc 
exTp{K X^A '''a), where the sum is over all the triangular pla- 
quettes (qubits) and ta = ±1 < when A £ E'. By inserting 
classical spin variables (T; = ±1 at the vertices and labeling 
the triangles A with triplets of vertices (ijk) we write P{E) 
as a partition function 



P{E) oc Z[K,t] 



(2) 
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TABLE I: Simulation parameters: L is the system size, Nss. is the 
number of disorder samples, tcq ~ 2*' is the number of equilibration 
sweeps, Tmin [Tmax] is the lowest [highest] temperature, and A^t the 
number of temperatures used. 
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12, 18 


20 


18 


2.200 


2.350 


31 


0.00 




24, 30 


20 


19 


2.200 


2.350 


31 


0.00 




36 


20 


20 


2.200 


2.350 


31 


0.02 




12, 18 


5 000 


18 


1.900 


2.400 


51 


0.02 




24, 30 


5 000 


19 


1.900 


2.400 


51 


0.02 




36 


5 000 


20 


1.900 


2.400 


51 


0.04 




12, 18 


5 000 


18 


1.700 


2.200 


51 


0.04 




24, 30 


5 000 


19 


1.700 


2.200 


51 


0.04 




36 


5 000 


20 


1.700 


2.200 


51 


0.06 




12, 18 


5 000 


18 


1.600 


2.100 


51 


0.06 




24, 30 


5 000 


19 


1.600 


2.100 


51 


0.06 




36 


5 000 


20 


1.600 


2.100 


51 


0.08 




12, 18 


5 000 


18 


1.400 


2.000 


61 


0.08 




24, 30 


5 000 


19 


1.400 


2.000 


61 


0.08 




36 


5 000 


20 


1.400 


2.000 


61 


0.10 — 


0.12 


12, 18 


5 000 


18 


0.750 


2.600 


38 


0.10 — 


0.12 


24, 30 


5 000 


19 


0.750 


2.600 


38 


0.10 — 


0.12 


36 


5 000 


20 


0.750 


2.600 


38 



Equation (|2]i is a 3-body classical Ising model with the cou- 
plings' sign given by t. When all ta = 1 the model is fer- 
romagnetically ordered at low T. Negative ta introduce frus- 
tration in the form of nets of domain walls. These can branch, 
a new feature not present in the random bond Ising model as- 
sociated with the TC. 

The relative importance of the different error homology 
classes P{E + Di) = Z[K,Ti] in ([T]i is given by the free 
energy cost of introducing a domain wall Di, because 



A,{t) ^ f3F{K,T,) ~ l3FiK,T) 



In 



Z[K,t] 
Z[K,T,] 



(3) 



The cost A; must be averaged over all coupling configura- 
tions, with p the probability for any triangle to have ta = — 1- 
Thus we are led to the study of a random 3-body Ising model. 
For low p and T (high K = J /T) the system is ordered 
and domain-wall fluctuations are suppressed: A; diverges 
with the system size for nontrivial domain walls. The criti- 
cal error threshold pc for error correction is recovered from 
the p-T phase diagram as the critical p along the N-line 

e-^-"^ =p/{l-p)M- 

Numerical details. — To determine the existence of 
a ferromagnetic phase we compute the finite-size cor- 
relation length il9tl . We start by determining the 
wave-vector-dependent susceptibility given by xi^) ~ 
(l/i^) 5^,^^.(S'j5j)TCxp[ik-(Rj-Rj)]. Here (• • • )t denotes 
a thermal average and the spatial location of the spins. The 
correlation length is given by 



U = (l/2)sin-i(g/2)v/[x(0)]av/[x('Z)]a 



1, (4) 



where q = {2n/ L, 0) is the smallest nonzero wave-vector and 
[• • • ]av represents an average over iVga disorder (error) sam- 
ples. Cm/^ ^ X{L^/''[T - Tc]), i.e., if there is a transition 



at T = Tc, data for ^m/i for different system sizes L cross 
at Tc [see, for example Fig.^a)]. The critical exponent f for 
the correlation length can be determined by a full scaling of 
the data, as shown in Fig. [3jb). We also probe the existence 
of a spin-glass phase by computing the spin-glass finite-size 
correlation length. 

The disorder in Eq. ^ increases the numerical complex- 
ity of the problem drastically with a behavior reminiscent of 
spin glasses ll20ll . To speed up the simulations, we use the 
exchange Monte Carlo method jlllj. Equilibration is tested 
by a logarithmic binning of the data. Once the last three bins 
agree within errors, we define the system to be equilibrated. 
Simulation parameters are shown in Table U 

Error threshold. — Figure [3] shows the temperature- 
dependent finite-size correlation length for different values of 
p. (a) Data forp — 0, the ferromagnetic case. The dashed line 
represents the transition temperature of the 2D Ising model 
Tc ~ 2.2692 I22II . The agreement with the numerical data 
is excellent, suggesting that coiTections to scaling are negli- 
gible, (b) Finite-size scaling analysis of the data in (a) using 
the exact exponent 1/ = 2/3. (c) — (h) Finite-size correlation 
length for different p values. For p ~ 0.108 marginal behav- 
ior appears and the determination of the transition is difficult. 
Because p = 0.107 shows a transition, and p = 0.109 shows 
marginal behavior, whereas p = 0.110 shows no sign of a 
transition, we conservatively estimate pc = 0.109(2) 1E3I1 . 
This is close to estimates for the TC where p'^'~^ has been con- 
tinuously improved from 0.1094(2) HI to 0.1093(2) iQ 
and 0.109187 11511 . The p-Tc phase diagram is shown in 
Fig. [11 the solid (blue) line being the N-line. We have also 
verified that there is no spin-glass order in the model (not 
shown). Finally, we ensure that our results do not violate 
the quantum Gilbert- Varshamov bound lfl3i[l4l 15, 24 ] where 
the encoding rate R{p) must satisfy R{p) < 1 — 2H{p), 
H(p) = —P log2 (p) — {1—p) log2 (1 — p) the Shannon entropy. 
For our estimate the bound is satisfied, since it 
lies under the zero-rate probability p ~ 0.110027. 

Conclusions. — In summary, we have computed the er- 
ror threshold for TCCs on a triangular lattice by mapping the 
problem onto a 3-body random Ising model on a triangular 
lattice. Using Monte Carlo simulations we find for the error 
threshold = 0.109(2) Q. Therefore, TCCs are as robust 
as the Kitaev toric code with the added benefit of being able 
to represent the whole Clifford group of quantum gates. The 
studied 3-body random Ising model highlights the relation- 
ship between spin-glass physics and information theory ll29ll . 
e.g., fully connected systems, and presents a new class of sys- 
tem exhibiting glassy behavior via 3-body interactions, with- 
out spin-reversal symmetry. Future work will focus on the im- 
pact of faulty measurements and the corresponding mapping 
to a (2 + 1) -dimensional random gauge model. 
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FIG. 3: (Color online) Finite-size correlation length (.m/L as a function of temperature T for different values of p. (a) p = 0. The data cross at 
the critical temperature of the 2D Ising model (dashed line), (b) Finite-size scaling analysis of the data for p = using — 2/3. The scaling 
is very good showing that corrections to scaling are negligible, (c) — (f) For p < Pc = 0.109 there is signature of a transition (data for different 
L cross) whereas forp > pc the transition vanishes [panels (g) — (h)]. 
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